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I. INTRODUCTION 


A. BACKGROUND 

Accurate ocean weather prediction is essential for successful fleet operations. The 
NANY needs superior data collection capability to obtain the data density and 
reliability necessary to produce consistently accurate forecasts and oceanographic data. 
However present prediction models are limited by the quality and quantitv of input 
data which come mainly from ships. The obvious approach to satisfy these purpose 
can be derived from satellite observations. Therefore, the NAVY planned the 
construction of the Navy Remote Ocean Sensing System ( N - ROSS ). [Ref. 1] This 
System consists of satellites which scan the earth surface and provide the fleet with 
timely worldwide knowledge of ocean data such as seasurface wind speed. wind 
direction, seasurface temperature, ice edge detection, ocean wave height and ocean 
photograpy. [Ref. 2] To satisfy the mission requirements the N - ROSS Satellite will 
carry several sensors. Among these sensors the Low Frequency Microwave Radiometer 
( LFMR ) is the most important and the most interesting from the dynamics of the 
spacecraft point of view. The function of the LF MR is scanning the earth surface and 
measure the seasurface temperature. To increuse the scanning area the deplovable 
reflector spins at [5 r.p.m. The sizes of this LF MR reflector and boom are relatively 
large compared to the N - ROSS Satellite itself. So the weight of this boom should be 
light, which makes the boom flexible. By this reason. there exist certain extent of 
deflection at the tip of the LF MR boom and this boom vibrates when this boom is 
spinning. Deflection and vibration due to this elastic deformation induce pointing 
error of the reflector in elevation and azimuth angle. However there is strict pointing 
Smor requirement of the LFMR. [Ref. 2] Therefore analysis of the flexible LFMR 


boom which supports the sensor payload is imperative for this research. 


B. STATEMENT OF PROBLEM 

The traditional approach to dynamics of a boom system is based on the 
assumption that the systems are composed of rigid bodies. Until recently, only a rigid 
body motion was assumed for the analysis. However, the flexible svstem includes a 
small elastic deformation as well as a large motion. These sinall elastic deformations 


include bending, twisting and axial extension. Development of a dynamic model 


1] 


including flexibility demands more accuracy for the svstem responses. Without 
considering these smail motion of a boom, we cannot expect a certain accuracy to 
maintain a spacecraft attitude and pointing control. (Ref. 3] Recently, efforts have 
been made to control maneuvers of mechanical systems which can not be adequately 
modeled using a rigid body assumption for all or some of the system components, 
especially in the fields of satellites, [Ref 4: pp. 257-264] 

Therefore, the development of a good dynamic model of a flexible system. an 
efficient dynamic equations formulation method and a good dynamic simulation schem 
are essential for the analysis of and identification of potential problems in the flexible 
LFMR system. 


C. THESIS OUTLINE 

In Chapter II, the development of an analytic model for the Lower Frequency 
Microwave Radiometer ( LFMR ) reflector boom in 3-dimensional motion 1s described. 
The large motion due to rotation is described by an equivalent rigid boom motion and 
elastic deflection of a flexible boom relative to the equivalent ngid boom motion is 
expressed using the mode superposition technique. The dyvnamucal equations for this 
model are formulated using the Lagrange’s method. 

In Chapter III, The computer implementations for the solution of the obtained 
equations are explained. For the modal analvsis of the system, NAsa STRuctural 
ANalvsis ( NASTRAN ) computer program was used and Dynamic Simulation 
Language ( DSL ) was applied to solve the simultaneous, nonlinear, ordinary 
differential equations. The LINPACK subroutines DGEFA and DGESL are also used 
in the dynamic simulation. 

In Chapter IV, simulation results are presented to investigate the deflection and 
pointing error of the LF MR. Comparisons are made bv changing the torque input 
condition. The problems considered are |) the effects of spin-up procedure on the 
pointing error of the LFMR reflector; 2) the effect of damping on the settling time of 
pointing error; 3) the equilibrium configuration of LF MR booms due to constant 
rotating speeds For the comparison purpose. two kinds of maternal. aluminum allov 
and composite material, were assumed as the LF MR boom materials. 

In chapter V, Conclusions are made from the research and some 
recommendations for future work in the area of the dynamic analysis of the flexible 


LFEMR reflector boom svstem are given. 


Il. FORMULATIONS OF DYNAMIC EQUATIONS 


A. INTRODUCTION 

In this Chapter, the dynamical equations of motion for the LFMR system which 
rotates in three dimensional space is developed. The 2-dimensional planar motion of 
the same boom was also studies and are presented in Appendix A. A simple dynanuc 
model of the LFMR system is developed for this analysis. The Lagrangian approach 
and the mode superposition technique are used for the formulation of the dynamic 


equations of the flexible LF MR system. 


B. DESCRIPTIONS OF THE MODEL 

The LFMR shown in Fig. 2.1 consists of four structures: a reflector, an upper 
reflector boom, a lower reflector boom and an electronics box. The reflector is attached 
to the top of the upper reflector boom. The upper reflector boom and the lower 
reflector boom 1s connected by a boom hinge. The electronics box is attached to the 
bottom end of the lower boom. 

For our analysis, the deplovable reflector is modeled as a concentrated mass at 
the tip of the upper reflector. We assume the boom hinge which connects the two 
booms is stiff and firm and there is no relative motion between the booms after the 
deplovment of the LF MR boom. Therefore we consider the whole system ( reflector, 
booms, boom hinge ) as a one body system. 
ite slr iiR svstem is connected to the Main Bus of the N - ROSS Satellite by a 
Spacecraft Boom. The attitude of the N - ROSS Satellite is controlled bv a Attitude 
Determination And Control Svstem ( ADACS ) verv accuratelv. [Ref: 5] Therefore. it 
is assumed that the spin axis and the base is reniain fixed in the reference frame fixed 
to the N - ROSS Satellite. The N - ROSS Spacecraft moves on a circular orbit with 
the spin axis always pointing the earth center. Hence, the gravitational force is in 
equilibrium: with the centrifugal force in the orbit plane: the LFMR svstem is in zero-g 
environmient. Therefore. in the dynamic model of the present studies, the reference 
frame fixed to the N - ROSS Spacecraft is assumed the Newtonian ( inertial ) reference 
frame and the LFMR system is in Zero-g environment. 

From the above assumptions the dvnanuc model of the LFMR svstem is defined 


as shown in Figure 2.2. The global coordinates X, Y, Z is fixed in the inertial reference 


ie) 


irame and a moving coordinate x, y, z ( local coordinates ), which is a body fixed 
coordinate system, !s defined as shown in Figure 2.2. The body fixed coordinate 
system ( local coordinates ) 1s attached to the base O. The local z-axes and the global 
Z-axes are the common axis of rotation. The local x, y-axes and the global X, Y-axes 


are in the same plane with angle difference 0. 


C. LAGRANGE’S EQUATION 

For anv system there must be same numbers of independent coordinates as the degrees 
of the freedom of the system to completely describe the motion of the system. The 
choice of coordinates is important in dynamic analysis. Such independent coordinates 
are called generalized coordinates and are denoted by the letter q.. For a system with a 
set of n independent generalized coordinates ie | 2. 3 See ,n), Lagrange's 


equations are expressed as [Ref. &7] 





d oT aie 6U 

eee eee = ee, (eqn 2.1) 
t  6q, éq, dq, 

(i = eee is vile) 


where T is the kinetic energy, U is the potential energy and Q, is the generalized forces 


which is defined as follows 


dR. 
6. = Va (eqn 2.2) 
© FOG; 


The dot over a variable means derivative with respect to time. F. is the force acting on 
particle ; and R; is the instantaneous position of particle 7 and may be expressed in 


terms of generalized coordinates 


R, = R, (Gy. Go, Vas «++ yl (eqniza) 


and OR; is the virtual displacement of the particle. 
To describe the motion of the LFMR boom system which composed of a large 
slow motion due to rotation and a small fast motion due to elastic vibration, two Kinds 


of generalized coordinates are defined. One is 9 for rotation of boom and others are q, 
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Figure 2.2 Modeling of LFMR boom system. 


eT n ) for h-th mode generalized displacement, where n is number of 
modes. 


Now Lagrange’s equation 2.1 is rewritten as 


d oT OT he! 


“ae ae) art ae 7 & —— 


and 


d goles oe OU 0 CRS) 
—[—]-—+—=Q eqinel, 
dt dq, dq, Cd, 

@ee— 1,2, 3,...-- ,n) 


Miecsextenmaiorce acting on the EF MIK system is assumed the torque t by a 
torque motor at the root of the boom. The contribution of this torque to the 
generalized forces is Q, = t and this torque does not contribute to Q,. Damping 
forces are assumed equal to the modal damping value in this analvsis. Since the modal 
damping values can be measured easilv. Therefore the contribution of damping forces 


to the generalized forces Q;, can be obtained using a dissipation function, 


I - e | 
) an > 26; @; M; q:*(t) (eqn 2.6) 
th 0 y 
with Fs ea 
h 6a, 


= — 26) @ My dh 


Where 
Sh | modal damping ratio of h-rh mode 
ne the natural frequency of A-zh mode 
My: the modal mass ot h-th mode 


Then equation 2.4 and 2.5 finally written as 


d | oT | oT ns ou 2.7) 

a — ——_ = 7t at 

dt — 08 08 o8 Jils 
and 

i OT ee . : 

at |, ) ~ “Gq, * Ja, ~~ 7 Fa eMa Sh le 


D. POSITION AND VELOCITY 

During the rotational motion of the LFMR system, the boom deforms. The 
deformed positions of a generic point in the system can be expressed the vector sum of 
a vector R(x) from the origin to the undeformed position of the point and a vector 


“= 


W(x.t) which 1s deflection, as shown in Figure 2.5. The notations show in the Figure 


} 


2.3 represent the following parameters: 


M : tip mass 

ae mass of electronics box 

u.: length of lower boom 

C, length of upper boom 

B: _ angle between two links 

ce applied torque 

Git: angular displacement 

A(t) : angular velocity 

R(X) : position vector of the point on the boom in the local x-direction 
ROR) position vector of the point on the boom after deformation 
W(x,t): deformation vector of boom 

ee unit vector of local x-direction 

ie unit vector of local y-direction 

Ke unit vector of local z-direction 

ig: unit vector of global X-direction 

ip unit vector of global Y-direction 


ky : unit vector of global Z-direction 


Then the position vector R(x.t) 1s expressed as 


R(x.t) = R(x) + W(x,t) (eqnezaa: 


a 


The undeformed position Rg(x)is represented by its components, 


R(x) = RX) + R(x) (eqn 2493 


gf 


= RA) ee Reco as 
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Figure 2.3 Parameters of the boom system. 


The elastic deflection W(x,t) is expressed as the modal sum as follows: 


W(x,t) = >) [ P(x) i + p(x) J + O7(x) k ] q(t) feqmn tel) 
I 
where 


1g 


@.*(x) iis 7h mode shape function in extension 


@."(x) j is ch mode shape function in translation 


From the equation 2.10 and equation 2.11, equation 2.9 can be rewritten in the form of 


R(x,t) = R(x) i + Rj(x) k (eqn Zoli 


+ VL ors) i + (x) § + 0700) k] g(t) 
1 


= [ R(x) + (x) g(t) Ji + TY Os) g(t) 11 
1 


+ [ Rx) + Y Ox) g(t) 1k 
1 


The velocity of the point in the Newtonian reference frame is obtained by taking the 
time derivative of equation 2.12 and applying the relation between the time derivative 


of unit vector i and j, 


1) 6 ae 
;=O0kK Xj= - Oi 


k=O0kxk=0 


the velocity 1s 


Rix.) = | = SO aaa ey ce q(t) i (eqn 2.13) 


+ [OR (xt OY (x) a(t) + ¥ (x) a(t) 1 


l 1 


+ LY ox) a(t) 1k 


) 


Eo NINETIC ENERGY, 

The kinetic energy of the svstem can be expressed by the summation of three 
different kinetic energies. One is the kinetic energy of the boom itself, another ts the 
Kinetic energy of the tip mass and the other 1s the kinetic energy of the R.F electronic 


box attached to the origin of the boom. 


20 


The kinetic energy of the bcom itself is 


1 fe : 
on { p ROG) e R(x,t) dm _ (eqn 2.14) 
llamo : 
la. J, R(x,t) ¢ R(x,t) pdx 
where 

dm: differential mass of the boom 

dx: differential length of the boom 

p: mass per unit length 

io: total length of the boom 

R(x,t) : absolute velocity of the points on the boom 


The kinetic energy of the tip mass ts 


I - : : 
i = > M R(£,t) = R(£,t) (edngiel) 
where 
M : tip mass 
ret) : the absolute velocity at the tip position 


Miae kinetic energy of the R.F electronic box is 





i = lepanee 16 
ea ? T57 (eqn 2 | ) 
where 
le mass moment of inertia for the electronic box 
oe2 2 ' 
OP: time derivative of angular displacement 


Thus the total kinetic energy of the svstem is 


Zi 


T= Tee il eee, (eqn 2.17) 


a ('R © R(x,t) d 
a , 0 (X,t) (X,t) mm 
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F. POTENTIAL ENERGY 

The potential energy U of the flexible boom system can be composed of the 
gravitational potential energy U, due to rotation of the system and the strain energy 
(or elastic energy) U, of the boom due to deformation. But in our model analvsis, we 
exclude gravitational acceleration and onlv consider the strain energy. The potential 
energy was determined by the work done by the static weight in the deflection. This 
Work is, of course , stored in the flexible flexible boom system as strain energy. 
In this thesis, We apply mide summation method to expand the deflection in terms of 
the normal modes of the svstem. The deflection of a boom without any external forces, 


satisfies the lollowing equation of motion. [Ref. 7] 


[ ELW (x,t) | piw(x,t) = 0 (eqn 20iay 
where 

De mass per unit length 

El: flexural rigidity 


and © represent the derivative with respect to x. The normal modes ;(x) of the boom 


must satisfv the equation 


[ Elp."(x) J) — 7 p p(x) = 0 (eqn 2.19) 


ae 


and its boundary conditions. The normal modes @;({x) also satisfy the orthogonality 


relation 


e — . 
J ,04(%) g(x) dm'= 0 (Gomes 1) (eqn 2.20) 


M: (for j=1) 
where M; is the generalized modal mass of the i-zh mode. 
As expressed in Appendix A, the deflection of the boom in the equation A.-+ for 


the general form 1s 


W(x,t)= )) (x) q:(t) (eqn 2.21) 


i 


and the generalized coordinate q(t) can be determined by applving Lagrange’s 
equation after setting up the kinetic and potential energies. 


Now the potential energy can be expressed as 
L= = ferw-2¢x t) dx (eqn 2.22) 
? 0 eke é 6 | a 
I eC ae af 
meee a, gs |) BIG. (x) @."(x) dx 
a: 0 J 


Multiplying the both sides of equation 2.19 by 0(X) and a{ter integration for the whole 


boom. equation 2.19 becomes 


) 


o>) 


C 
J O)(%) LEL 9;(x) | dx (eqn 2.2 


>,¢ 
= OD: J oil) @}(X) dm 
After integration by parts and using the boundary conditions, equation 2.23 becomes 


e 
j oF! [ @;°(X) OF (x) J dx (eqn 2.24) 
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( 
=e oa 
= 07 J (x) O(x) dm 


From the orthogonaiity condition, the equation 2.20 and from equation 2.24, the strain 


energy in equation 2.22 becomes 


U = ¥ 7 M; q(t) (eqn 2.25) 


G. DERIVATIONS OF EQUATIONS OF MOTION 
From the equation 2.13, the dot product of R(x,t) 1s 


R(x.t)eR(x.t) = [ — OY g(x) a(t) + Y 92(x) a(t) F (eqn 2.26) - 
l l 
~[O9R(WWt OV oq + Vey 40 7 
1 ; 1 
Wea 2 
+ EL 9s) a(t) ]° 


] 


= 0215 ox) Qo + 1D 989 Qo? 
i 1 
26 


¥ (x) a(t) © () q(t) 
1 l 


+ R(x? 92 a 92 = 9.*(Xx) g,(t) \- | y 2 (x) a(t) \2 
i 


+ 2R(x) 07 q(t) + 28 F Ox) a(t) F 2x) qt) 
1 


+ IROL G(x) a(t) + 1E 900) a(t 7 
l 


= 1 { Yes) a(t}? + (Lew a(t}? + (¥ e200 g(t) 37] 


1 1 1 


an 


$) , P; ; =) 
+ [EX Ox) g(t) }° + CN oPr(x) a(t) Po + CX 0200) g(t) Jo] 
l l 


1 


aa? X(x) g(t) F @X(x) a(t) + RL2(x) 92 
1 


ee ex) 6 » Ox) aft) + 2 ; ~ .*(x) q,(t) ~ 0.¥(x) q.(t) 
1 
+2R, 0) 9Y oF Nix) q(t) — 02 (E err 


1 


Substituting equation 2.26 into the equation 2.17 and apply orthogonal relationship 


e | a | 
) ot PEER) OO) + FOO GFE) + O70) 70) } dm (eqn 2.27 


+M (OD OH +O7(H OH + OH OE) } = 0 (for ij) 


Mi; (for 1=j) 
then the kinetic energy is 


e 5 l 
T = — g° V g.7(t) M, + ae g(t)” M. (eqn 2.28) 
1 


; 5 ul 
— OD Lait) (OF) O79) Ox) dm + M OF(6 GAO | 
i 


l 


») 
ao 





° e 
2 2 
ae |i i R,-(x) dm + M R,(8) 


ic 
+9-Vv a(t) { J R.(x) @.*(x) dm + MR,(€) ,*(0) 
=A oF 0 xX 1 Xx 1 


° = a ié . 
+ OYE at af OM) 07x) dm + M 9,%(6) O70) | 
| 


: 4 
+ OLY a) (J Rs) 9260) dm + MR) (0) | 
1 
- 8 i CE 970) a(t) 2 + MEY O48) a(t) }71 
1 1 
+—1, 0 
D977 


aS 


By applying the kinetic energy and potential energy expressions ( ean 2.25 ) and ( eqn 
2.28 ), to the Lagrange’s equations 2.7 and 2.8 then equations of motion will be 


reduced as follows. The detailed derivation process can be formed in Appendix B. 


B10 M, ae R, 2¢x) dm + M R, 2(e) (eqn 2.29) 


7 
+ 2N att) CJ ROY O09) dm + MRO O;(0 } 
1 


4 


— 2{§ (YS oF@ al)? + MCE Ob a(t)? + 1 | 
() 1 1 1 a7. 

+ 26F diy [ ait M+ J RCO 980) Um + MRD 9. (£) } 
1 

=3 


E ato 4, OX) Ox) dm - M 0%) p70) } ] 
j 


SEO on 
l 


if 
i (p."(X) p, (x) dm + M 9,*(€) 0; *(€) 


e 
(p.”(X) @, (x) dm — M 9,°(0) 0, (Ga 


= ‘ 
- oR (x)dm — M R,(€) 99(€) ] = 
Jt 
qy(t) Mee 0(y q.(t) Ou p.*(x) @ (x) dm + M @,*(€) 0, 9(€) (eqn 2.30) 
1 


c 

= { ; p.” ( x) 0p, ~ydmi = vi p.” (C) (8) } 
"4 

- J RO p(x) dm + M R,(€) @,7(8 | 


+20¥ WOLL X(x) @9(x) dm + M @,*(€) 9,*(6) 


— 


1 


e f 
— J,o209 @, 809 dm — M @5(8) O50] 
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4 e 
= 7 (gi) M; — J RCO) M(x) dm — MR,(E) ONE) 


t 
+ 22 a(t) Cf 0,700 9,59 dm + M O12) 9,78) 3] 
1 


+ 0°, M, q(t) = 0 


Now let’s define the following quantities: 


me 
Mg = [Lat M, + J R(x) dm + M R7(0) 
1 


e 
+ 2D a(t) (J RC) OM) dm + MRO) OO } 
1 


oul 
-26f (Loxwaor?+ MC Lea? +t 
l 


| cmediiel 


a 
Méaq. 2 [q.(tjae, + J Rx®) o.*(x) dm + M R,(€) @.*(@) } 
oe 
— 2¥ a(t) J 9,) 0708) dm - M 0,76) 018 } | 
| e 
Mg = ~[LaH (foo) oF(x) dm + MoO o%(6) 
i 


e 
— 07) @*(x) dm — M 92(€) @,*(€) } 
e : : 
— J R(x) 92(x) dm — M R,(€) (8) J 
e | 
Mon = 0X at) CJ OMX) O00) dm + M O(8) (0 
1 


4 
Bs 9209 @,,*(x) dm - M @,%(£) @, *(€) } 


4 
2 f Rw 0, 9(x) dm + M R,(£) 9, *(€) | 


a 


4 
Méq, = 2[] Or) Q(x) dm + M .%(€) @,7(0) 


£.. 
= J &) O, (x) ém — M 9."(€) Oy (8) |] 


e 
Fo = O° Lait) M, — J RAx) 9x) dm — MRO) O;%(0) 


4 
+ 2 » q.(t) { J e700) @, (x) dm + M @*(€) 6,°(8) 3] 
i 


aes = M O,7 q,(t) 


Then equation 2.29 becomes 
ee « < a . | eat _ ; 
Mg 8 + MQ q, 0) q(t) - My St oy (eqn 2:ciy 
and quation 2.50 becomes 


Mg, 8 + Mea, O0Pq (t) erie (tee Fo ge Cees (eqn 2.32) 
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Il. COMPUTER IMPLEMENTATION 


A. INTRODUCTION 

To solve the equations of motion, two computer program, Dynamic Simulation 

Language ( DSL ) and NAsa STRuctural ANalvsis ( NASTRAN ) program are used. 
NASTRAN is a general purpose digital computer program for the analysis of large 
complex structures. [Ref. 8] This finite element computer program was used in the 
modal analysis which determine the mode shapes, generalized modal mass, generalized 
stiffness, natural frequency of the LFMR model. Then these properties directly 
inputed to the DSL program to get a set of solutions. 
DSL is an IBM/VS FORTRAN-based simulation language for digital simulation of 
continuous system. [Ref. 9] It is one of the most effective for the solution of 
continuous modeling and simulation problem with computational power (automatic 
double precision and accurate timing), whether the problem is time based or not. 

For the integration method to solve simultaneous nonlinear second order coupled 
ordinary differential equations. Runge-Kutta method was chosen. Runge-Kutta fifth 
order integration method ( RKS5S } is self-starting, stable and automatically determine 
the step-size but this needs excessive computer time. 

Two different cases are analyzed. One is the LFMR boom made of Aluminum 
Alloy and the other is the LFMR boom made of Isotropic Graphite Epoxy composite 
material ( 1300/5028 ( 0/90'45/-45)¢ ). The mass distribution of the two cases are the 
same. Therefore the boom model of graphite epoxy composite material was stiffer since 


the mass density was linear. 


B. MODAL ANALYSIS 

As previously mentioned, we consider our model as a one body system and we 
dont need any compatibility conditions which is necessarv in multibodv dynamics. For 
the modal analvsis we equally divide the whole boom with fourteen grid points. Figure 
3.1 shows Finite element model of LFMR boom. As the number of grid points 
increases We can get more accurate results. But Degrees of Freedom ( DOF ) of the 
System also increases. Fourteen grid point is sufficient for our our analvsis. 

In the modal analvsis of the boom . Modified Given’s method ( MGIV ) was 


applied and for the purpose of simplifving the equation of motion we normalized the 


78) 


moce shape such that the generalized modal mass equal to unity. From: the 
relationship of the generalized stiffness, generalized modal mass and natural frequency, 


the generalized stifiness K; is 


2 


and reduces to 


a 


Kj = OF 


Only the first 2 modes were needed in the 3-dimensional dynamic analysis with 
sufficient accuracy. 

NASTRAN program is shown in Appendix C and some outputs are tabulated in 
Table | and 2 for 3-dimensional motion. Figures 3.2 and 3.3 show the first and the 
second mode shapes for three dimensional motion. In Figures 3.2 and 3.3, the left side 
figure is the projection to x-z plane and right side figure is projection to y-z plane. 

The first mode shape shows there is no in-plane vibration and the second mode 


shape shows there is no out-of-plane vibration. 


Tables 
REAL EIGENVALUES OF ALL MUANUIOAELO abs 


mode eons ae generalized generalized 
reKey mass( MI.) stiffness(K, ) 


Jee aI5E +00 | 53866236 -01 586623E re / 1.0000008+00_| lf p 114sa9ge tor 


1 .QQOOOOE + 00 1.273648E+01 






D.00ss2 (sa 00 





C. DSL PROGRAMING 

DSL was implemented for the solution of a set of simultaneous, nonlinear, 
second order, ordinarv differential equations. DSL offers a new simulation tool that 
speeds problem analvsis and solution for a wide variety of application. It is adopted for 
simulation because it requires less skillful time for problem analysis and provides 
quicker, more comprehensive plotting through GRAFAEL and sophiscated line or 
print plot. These advantages can translate directly into higher productivity and cost 


Savings. 
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TABLE 2 
REAL EIGENVALUES OF COMPOSITE MATERIAL (3D) 


mode i oe Fa generalized 
nO, mass( stiffness(K, ) 


Bi ja | 4.287485E+00 OO | 6.823744E-O1 §23744E — 01 1 .QOOQOOQVE + 00 Sole! 









ve — 1.000000E+00 | 1.990350E+01 


DSL is a high-level continuous simulation language which incorporates VS 
FORTRAN as a subset. Because of its tremendous flexibility, DSL facilitates the 
solution of nearly all priblems involving time-dependent differential equations. Thus 
DSL readily assists in the dynamic analysis of transient behavior of dvnanuc systems. 
Also DSL 1s easily learned and applied to many problems in science. engineering, 
mathematics and management. Coding is simple. execution is rapid and results can be 
displayed graphically. For any one involved in simulation modeling DSL offers 
increased power for faster problem solving and the user choice of nine integration 
method ( fixed-step, variable-step. variable-step & variable-order method ). 

In our miodel analysis, Runge-Kutta fifth order method ( RKS ) was used because 
it is self-starting, stable and provides good accuracy. To code the equations of motion. 


the equation 2.31 and 2.32 have to be rewritten as follows: 


Mg 8+ Mg q(t) = t — Mg, 9¥ av (eqn 3.1) 


and equation 2.30 becomes 


Mo, 9 + My ait) = Fo — Fe — Mga, 9 ¥ 4, 


1 


(t) (Commo ) 
(Clas Be n ) 

Two ways solving these systems of equations are shown. One uses matrix algebra and 

utilizes subroutines from the library LINPACK, the other uses subroutine REGULRA 


provided by DSL. For our system analvsis we select matrix algebra because it is easv 


to use for large numbers of variables. 


31 





Figure 35.1 Finite element model of LFMR boom. 


a2 





Figure 3.2 First mode shape. 
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Figure 3.3. Second modesinaee 


Using matrix notation, the equations 3.1 and 3.2 may be rewritten as 
[M]{X} = (f} 


where the matrices M, X, and f are defined as follows: 


Mg Mg ---- sere Mg 
lo, 0 
lee Doo soe 0 
| \10 i 
if peor, 


Vlg, CMIARE Sa rce  er ay Gur car a aaa era M 
) 
qy 
- 2 
‘Xi = 
a | 
el J 
Mgq; 
Fe — Fe, — Mogq, 9 qy 
(f} = 


Then 
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(X}=(M]7!(83 


and by successive integration of { X }, the solution vector may be obiained. Since the 
elements of the [ M ] and { f } matrices are time dependent, the { x } vector must be 
computed at each integration step. A matrix decomposition subroutine ( DGEFA ) 
using Gaussian elimination and a subroutine which uses this decomposition to solve a 
matrix equation ( DGESL ) are called directly from the LINPACK library) Becarae 
DGESL returns the solution vector { Xx } in the right-hand-side vector f, the vector 
name { X } does not appear explicitly in the program. 

The coefficients of [ M ] and { f } matrices are calculated in each time step. 9, qi 
(ic dee. eee , n ), deflection and slope at the tip position in each direction were 
calculated at each time interval of interest. The graphs of these variables vs. time were 


also obtained. The computer program for dynamic analysis is coded in Appendix D. 
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IV. RESULTS AND DISCUSSIONS 


Computer simulation was done in four areas for double link boom system in 
three dimensional motion to investigate the equilibrium configuration and the vibration 
amplitude of the tip position. The effects of the spin up procedure on the pointing 
error of the reflector is studied: the first analysis was the comparison of the three 
different torque histories by maintaining magnitude of torque until rotating speed of 
boom reaches to 15 r.p.m. Secondly the magnitude of applhied torque was changed in 
three cases for one of three torque applying methods above. The effects of structural 
damping of the boom on the settling time is also studied by changing the modal 
damping values. The magnitude of the deflection and slope at the equilibrium 
configuration at three different rotating speeds are investigated by simulating free 
motions of the system with the three different initial rotating speed and undeformed 
configuration. Also some comparisons were made for double link flexible boom svstem 


in planar motion in Appendix B. 


A. MATERIAL PROPERTIES AND PARAMETERS 

For the comparison of the computer simulation results, two kinds of material were 
chosen. One is Aluminum allov 6061-T6 (99% Al-1% Mg) and the other is an 
Isotropic Graphite-Epoxy Composite material (1300/5208 (0/90/45/-45)¢). [Ref. LO] 
Same mass, same length of each link, same outer dianzeter was used for these two 
fiiaeenials. Size of the electronic box is | ft. cube. The properties of the materials and 


some geometric parameters are given in Table 3 and Table 4 respectivly. 


fee EFFECTS OF TORQUE APPLYING PROCEDURE 

Three cases of torque applving method were chosen without damping for 
comparison. To maintain 15 r.p.m of rotating speed, we select the maximum magnitude 
of torque is 10 lbs-in. This value may be appropriate to maintain 15 r.p.m in 20 to 30 
seconds. 

In the first case, 10 lbs-in torque was abruptly applied from the begining. This 
method is unlikely in the real situation but we chose this case for comparison purpose. 
In the second case, torque was parabolically applied until reaching 15 r.p.m then cut 
off abruptly to 9. Finally torque was applied same manner as case two but linearly 


decreased to 0. Figure 4.1 shows three different spin-up schemes. 


al 


TNBIGES 
BOOM PROPERWIES 


PROPERTY 


ALUNU NU SE NEE Om 


COMPOSITE 





specific weight ( 7 ) 0.0980 Ib‘in 0.0588 lb in? 


modulus of elasticity ( E ) lO. TE = 0Gtns; lO.TE + OGaes 

modulus of rigidity ( G ) S1E+ 06S 
E 
Sr 


thickness (t) 0.1344 in 0.2319 in 


cross area ( A ) 1.21017 2.0168 ine 


| mass per unit length ( p ) 3.0690E-04 Ib. in 3.0690E-04 Ib. in 
jane) 


polar moment of inertia ( J ) 2.4894 in? 3.8902 int 


In these three cases angular displacement 0 was changed parabolically as shown 


250 
3.0 in 
2.5362 in 





in Figure 4.2 while torque is aprlied, then linearly increased as torque was cut off. The 
cases B and C are almost indentical and overlapped in the [igure 4.2. Angular 
displacement was equally increased with time for aluminum alloy and composite 
material because total svstem mass is same. As shown in Figure 4.3, angular velocity @ 
varied linearly with some oscillation in the begining for the first case because torque 
was suddenly applied, then it gradually decreased until reaching 15 r.p.m. After cutting 
off applied torque, 0 oscillates periodically because the boom deflects to the positive x- 
direction and the concentrated miass center moves to the rotating axes, this means 
radius of rotation decreased and moment of inertia also decreased. From the 


conservation of angular momentum, 


Oe 


Ho =mr = constant 
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TABLE 4 
GEOMETRIC PARAMETERS 












PARAMETER VALUE 


Posse [ee 


; tip mass 37.5 Ibs 


So as the tip position changes, 0 varies with reciprocal of r~. Therefore angular 






C, ; length of lower boom 








velocity oscillates periodically. 
. 

In case two and three, 8 increased smoothly without oscillation until reaching 15 
r.p.m then it had relatively small oscillation as soon as torque was removed gradually. 
Figure 4.5 shows the angular velocity of the Aluminum Alloy boom in magnified seale. 
For the composite material as shown in lrigure 4.4, magnitude of oscillation was much 
smaller than that of aluminum allov so it looks no escillatron. But magnified angular 
velocity of constant part shows obvious oscillation in ligure 4.6. ]:.ures 4.7 - 4.10 
show the variations of the generalized coordinates q, and q, respectively. 

During rotation deflection in x-direction was dominated and deflection center 
reaches to its equilibrium position in x, Z-direction then oscillates harmontieally as 
shown in Figures 4.11 and 4.12 but displacement in y-direction was rapidly inereased at 
the begining then decreased gradually to 0 and as soon as applied torgue was removed 
it oscillates harmonically as shown in Figures 4.13 and Figure 4.14. 

Characteristics of displacement and slope at the tip position its verv similar to its 
generalized coordinates, and the generalized coordinates q, in Figures 4.9 and 4.10 was 
dominated for the system deflection. Displacement in z-direction was shown in Figure 
4.15 and 16 for two material booms. Displacement and slope change as shown in 


figures 4.17 - 4.22 was apparently different from case A and case B but alter cutting 


39 


off the torque its characteristics of deformation was not much different. For the case 
three deflection and slope increases smoothly while torque is applving then as torque 
was removed it reaches its equilibrium position and oscillates with small fluctuation 
ero 

Figures 4.11 -4.24 show displacement and slope at the tip position in each 
direction for aluminum alloy boom and composite material boom. 

As We mentioned previously, we are mainlv interested in the pointing error at the 
tip position in both elevation and azimuth. For the composite material boom, all 3 
cases Were satistied the pointing error requirement ( elevation and azimuth angle error: 
0.0328° ). But for the Aluminum Alloy boom, torque applying case A and case B were 
not satisfied the elevation angle error requirement. Figures 4.19-4.22 show the 
characteristics of elevation and azimuth angle variation, and in Figures 4.23 and 4.24, 
the effects of torque applying procedure are compared. From these results, pointing 
error at the equilibrium position depends on torque removing procedure. As we see in 
Figures 4.23 and 4.24, pointing error in case A and case B is almost same but it is verv 
different in case C and pointing error varies more sensitively with flexible material as 
torque removing procedure varies. Consequently, if we applv the torque more gently 
and remove slowly with sufficient time, the pointing error both deflection and the slope 


can be reduced to much smaller values. 


C. EFFECTS OF CHANGING MAGNITUDE OF TORQUE 

In this section we investigate the effects of the maximum magnitude of applied 
torque. For the simple comparison, we select second case of torque applying method 
used in previous section. Magnitudes of torque chosen are [0 lbs-in. 20 Ibs-in and 40 
lbs-in as shown in Figure 4.25. As the magnitude of applied torque increases 4 times. 
applving time was reduced to almost one third. All characteristics are same as the 
second case of previous section for 3 cases shown in Figures 4.26-4.36 but magnitude 
of pointing error increased linearly as the maximum torque increases as shown in 
Figures 4.37 and 4.38. Figures 4.26-4.43 show angular displacement. angular velocity, 
generalized coordinates, displacement and pointing error for the boom made of 


Aluminum Alloy. 


D. EFFECTS OF DAMPING COBLRICI IES 


In this section, we will investigate the damping effects for settling down the 


oe 


vibration of the boom. As we mentioned before, in the 3 cases of torque applying 


40 


procedure, the second case did not satisfv the pointing error requirement in elevation 
angle. We will not mention about the first case because it is unlikely to be used in real 
situations. For the comparison purpose, we arbitrary choose 3 modal damping 
coefficient such as 0.2 %, 0.5 % and 1.0 %. To see the results more clearly, we 
simulate the model in 200 seconds. Figure 4.39 shows the spin-up procedure used in 
this analvsis. Figure 4.40 shows vibration settling down ratio during the first 200 
seconds. In this Figure, it is hard to recognize the effects of all 3 cases but in Figure 
4.41, we can clearly see the vibration settling down for each case. 

Two trials were made for investigation. Firstly, we tried to find how much 
vibration was settled down in 200 seconds for each damping coefficient. Secondly, we 
found how much time was required to meet the pointing error requirement for each 
case. Figure 4.42 shows elevation angle pointing error decreasing ratio in 200 seconds. 
Each value stands for magnitude of pointing error ratio to its initial value which is the 
magnitude of pointing error without damping. As shown in Figure 4.42 pointing error 
decreased exponentially as damping coefficient increases linearly. Figure 4.43 shows 
desired time to meet pointing error requirement for each damping coefficient. As 
damping coefficient increases linearly, desired time decreased exponentially. 

THE EQUILIBRIUM CONFIGURATION AT DIFFERENT ROTATING 

SPEEDS 

Initial conditions were given for 3 cases such as 5 r.p.m, 10 r.p.m and 15 r.p.m 
with damping coefficient of 2 %. From this result, we could investigate the magnitude 
of deflection in each direction and slope at the tip position. Figure 4.45 shows angular 
velocity ( rotating speed ) variation with time. As the initial rotating speed increases. 
magnitude of oscillation was also increased and the rotating speed at the equilibrium 
state also increased. As we mentioned in section B, angular momentum is constant in 
the system since no external torque is applying. Therefore as the boom deflects to the 
positive x-direction, radius of rotation decreases so angular velocity increases 
consequently. Figure 4.45 shows elevation angle change and its magnitude of 
fluctuation change. As the magnitude of initial speed increases linearlv, elevation angle 
change and magnitude of fluctuation increases exponentially. This result is shown in 
Figure 4.47. Magnitude of deflection in x-direction and z-direction was also increased 
as the rotating speed increased as shown in Figures 4.47 and 4.48 respectively. This 
effects was more sensitive for the flexible material boom. The oscillation of azimuth 


angle centered at 0. 
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Figure 4.1 


Applied torque vs. time. 
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Figure 4.2 Angular displacement vs. ume. 
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Figure 4.3. Angular velocitv vs. time ( AL. ). 
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Figure 4.4 Angular velocity vs. ume ( COM. ). 
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Figure 4.5 Magnified angular velocity vs. time ( AL. ). 
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Figure 4.6 Magnified angular velocity vs. time ( COM. ). 
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Figure 4.7 First mode generalized displacement vs. time ( AL. ). 
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Figure 4.8 First mode generalized displacement vs. time ( COM. ). 
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Figure 4.10 Second mode generalized displacement vs. time ( COM. ). 
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Figure 4.11 Displacement in x-direction vs. time ( AL. ). 
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Displacement in x-direction vs. time ( COM. ). 
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Figure 4.13. Displacement in y-direction vs. time ( AL. ). 
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Figure 4.14 Displacement in v-direction vs. time ( COM. ). 
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Displacement in z-direction vs. time ( AL. ). 
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Figure 4.16 Displacement in z-direction vs. time ( COM. ). 
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Figure 4.17 Mlagnitude of deflection at tip position vs. time ( AL. ). 
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Figure 4.18 Magnitude of deflection at tip position vs. time ( COM. ). 
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Figure 4.19 Elevation angle change vs. time ( AL. ). 
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Figure 4.20 Elevation angle change vs. time ( COM. ). 


6] 


ee 
Tn ee ee ee 

SS —_——re 
eee 


an Sa 


oa, 


A 






5 ie IE 
oe ERR sen 





Ee 
,_ ——_—_—— 


pe — 
SLE 
ue 
: SS O_—H>O>-.”—__S > 
< —— > 


eee 
, pS 


vy - 


*. 
= SF e t - 


I3T-OT SINGH J19NG HINWI ZY 


Figure 4.21 Azimuth angle change vs. time ( AL. ie 
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Figure 4.22 Azimuth angle change vs. time ( COM. ). 
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23 Pointing error in elevation angle vs. torque applying procedure. 


Figure 4. 
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Figure 4.24 Pointing error in azimuth angle vs. torque applying procedure. 
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Figure 4.25 
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Applied torque with changing magnitude vs. time ( AL. ). 
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Angular displacement vs. time ( AL. ). 
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Figure 4.27 Angular velocity vs. time ( AL. ). 
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Figure 4.28 
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Magnified angular velocity vs. time ( AL. ). 
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Figure 4.29 First mode generalized displacement vs. time ( AL. ). 
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Figure 4.30 
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Second mode generalized displacement vs. time ( AL. ). 
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Displacement in x-direction vs. time ( AL. ). 
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Figure 4.32 Displacement in y-direction vs. time ( AL. ). 
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Figure 4.33 Displacement in 2-direction vs. time ( AL. ). 
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Figure 4.34 Nlagnitude of deflection at tip position vs. time ( AL. ). 
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Figure 4.35 Elevation angle change vs. time ( AL. ). 
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Figure 4.36 Azimuth angle change vs. ume ( AL. ). 
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Figure 4.37 Pointing error change in elevation angle vs. magnitude of torque ( AL. ). 
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Figure 4.59 Applied torque vs. time with damping ( AL. ). 
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Figure 4.40 Elevation angie change vs. time with damping ( AL. ). 
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Figure 4.42 Elevation angle error decreasing ratio in 200 sec. vs. damping coefficient. 
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Figure 4.45 


84 


(93S) Wil 





001 08 09 O» 02 0 
1°0 
0 
“0 
70 
8°O 
. a 
ga a a ae eee NTN NIC NR |B 
| 
— —— — Wl SI 
ee Wedd Ol $ 
Hdd S 
AKARNAAKAN S| \ Wr" 
NN GRU of 1] NG aN Na, \ 
a alk vy SMM 


“ha 


( 338/084 } ALIDOIZA sd ONY 


Initial angular velocities vs. time with damping ( AL. ). 
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85 


Ni ——— -—— 


15 R.PLM — — — 


So Rae 
DOR at: 


MME? == Se eG ee, ee EE 





60 60 100 


TIME (SEC). 


40 


( 
| 

| 

g 
2 
S 
° 
2 
S 
SS 
<= 


20 


€ 930 ) SONYHD STSNY NO! LYAR13 


Figure 4.45 Elevation angle change with damping ( § = 2 % ) vs. rotating speed ( AL. ). 
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Figure 4.46 Elevation angle change at the equilibrium position vs. rotating speed ( AL. ). 
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Figure 4.47 Deflection in x-direction vs. rotating speed ( Ale: 
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Figure 4.48 Deflection in z-direction vs. rotating speed ( AL. ). 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The purposes of this research are to investigate the the effects of the flexibility of 
the LFMR boom on pointing error of the reflector in elevation and azimuth angle and 
to identify the parameters which are important in pointing error. In this studv, the 
effectiveness of the application of Lagrange’s method and mode superposition 
techniques, computer simulation techniques are also investigated. 

The results indicated that the Lagrange’s method and mode superposition 
technique was very effective for this double link boom model analvsis. The effects of 
flexibility of boom was sufficiently expressed with few mode and the boom flexibility 
Was verv affective in pointing error. As the material becomes stiffer with light weight, 
its pointing error becomes smaller. The torque applving and removing procedure 1s 
very important to the magnitude of pointing error. [f we apply the torque graduailv 
until reaching desired rotating speed and remove slowly with sufficient time, we can 
reduce the pointing error to much smaller values within requirement. The vibration 
Settling time decreases exponentially as modal damping coefficient increases and the 
pointing error is linearly dependent on the magnitude of applied torque. The deflection 
in each direction and elevation angle change in equilibrium condition increase 
exponentially as the rotating speed increases. Sensitivity of deflection depends on 
flexibilitv of the boom. 


B. RECOMMENDATIONS 
In this researcn, we regarded the deplovable reflector as a concentrated tip mass. 
In this case, we dont need consideration of the flexibility of the reflector but in actual 
model the flexibility of the reflector will affect the accuracy of the pointing error, 
therefore flexibility of the deplovable reflector has to be considered for the future work. 
As we have seen in the comparison of two material boom, the flexibility of boom 
was very affective in the pointing error. Stiffer material with light weight will play 


important role in the reduction of pointing error of the flexible boom. 
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APPENDIX A 
DERIVATIONS OF THE EQUATIONS OF MOTION 


In this Appendix, we develop simplified dynamic equations of motion for 2 cases 
of the boom system with tip mass. 

The large motion caused by rotation and the small motion created by elastic 
deformation will be expressed bv generalized coordinates then applv Lagrange’s 
equations to develop the equations of motion. I[n these analysis, the local rotarv 
inertia and shear deformation of booms are negiected. 


From the equations 2./ and 2.8 Lagrange’s equations are 


d | OT OT r. G \1) 
-— =f eqn A. 
a 60) 8tCD “aa 
and 
d | OT ) OT OU , AD) 
——|——- | - — —— = Can sr.e 
dt dq, dq, 44, 
Mme 1, 2,3,..... <n) 
where 
Q: the generalized coordinates of the large motion 
qn: the generalized coordinates of the small motion 
oe applied torque to the system 
ig number of modes (or number of degrees of freedom) 


I. SINGLE LINK BOOM SYSTEM IN PLANAR MOTION 
a. Geometry of the system 
In model analysis, since the extension deformation is negligible, only the 


bending deformation is considered in the analysis. 
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Figure A.1 Parameters of the single link boom system. 


tip mass 

mass of electronics box 

length of the link 

applied torque 

angular displacement 

angular velocitv 

position vector in local x-direction 

deformation vector of boom 

position vector of the point on the boom after deformation 
unit vector of locai x-direction 


unit vector of local y-direction 


x { 
xq 


eK: unit vector of local z-direction 
unit vector of global X-direction 
ip: unit vector of global Y-direction 


Ko : unit vector of global Z-direction 


b. Position and velocity 
During rotation the boom deforms and the positions of any point in the 
system can be expressed the vector sum of a vector R,(x) from the origin to x and a 


vector W(x,t) caused by deformation. Then the position vector R(x,t) is expressed as 
R(x,t) = R(X) + W(x,t) (eqn A.3) 


where R,(x) is only x dependent variable 


SO 
R(x) = R(x) i 


and W(x.t) is deformation obtained from modal summation method 
then 


W(x.t) = > o,(x) q,(t) (eqn 2.21) 


1 


where 


(x) is i-th mode shape function 


q(t) is -7h mode generalized coordinates 


If we consider the deformation of the boom consists of translation and extension. then 


equation 2.21 becomes 


Wix,t) = > [ O%(x) q(t) i + 92(x) 4.(t) §] (eqn A.4) 


1 


where 
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p(x) Lis 2h mode shape function in extension 


.*(x) j is zh mode shape function in translation 


From equation 2.21 and A.4 the position of the point on the boom is 


R(x.t) = R(x) i + VY [ @X(x) i + @F(x) j] a,(t) 


i 


(eqn A.5) 


Now the velocities of the point on the boom is necessary to formulate the 


kinetic energy to apply Lagrange’s equation. The velocity is obtained bv simply 


differentiating the equation A.5 
then 


Rixt) = Rix)i + VY (a(x i + O89) $1 g(t) 


1 
+Tleswit of ila(d 
i 
={R, + Y ox) a(t) ]i + Y O20 a | 
1 i 
= NOS) g(t) i + OP (x)a( GF 


i 


But time derivative of unit vector i and j are 


=Ok x i= 6 
j= Ok x j= — Oj 
k=O0kxk=0 


substitute these quantities to the equation A.6 and simplify 
then 


R(x,t) = {R. + & (x) q(t) ) 07 — O Oras 
1 i 
+ ¥ p(x) g(t) i+ Vps(x) g(t) j 
1 1 
= { ¥ (x) q(t) — 6 © @3(x) q(t) ] i 
1 1 
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(eqn A.6) 


(eqn A.7) 


+ [| R(x) 8 + 6 ¥ ox) q(t) + ¥ p(x) g(t) ]j 


c. kinetic energy and potential energy 


Recall the equation A.7 and drop the all terms concerning .*(x) and express 
p."(x) simply @.(x) 


then equation A.7 reduces to 


R( x.t) 


ll 
ye 
De 
om < 


— OY (x) g(t) i + SY 0200 a(t) Jj (eqn A.3) 
i i 


SeomGaiax)g(t) i= | R(x) 8 2 (x) a(t) 11 


From the equation 2.17, the kmetic energy of the whole system is 


T= ae oa Meee - Te eon) 
= C : 2 
2 i R(x,t) ¢ R(x,t) dm 
l ° ae 
3 > Mi R(€,t) 4 R(£,t) 
” 
but 
R(x.t)eR(x.t) =[— OT 9 (x) q,(t) } +R (X) a S (x) q.(t) }° (eqn A.10) 
i i 
= PLY a(x) ao F + [Roo 8)? 
I 
+ 2R(x) OY (x) a(t) + LY 9x) Gt) 
1 
and 
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IE as) ao “24 (x) 70) + 2 E 0x) Ox) a,(0a,(t) 
i 


2 g(x) §(t) 2 = F aft 820 + 25 (0) 1) G00 


(for i#j;2=0) 


Substituting into equation A.10 


then 


R(xeR(x.t) = (Le 70 a7) + 2D) OW) ada(t) — (eqn ALI) 
i ij 


(x) 87] + 2R Gx) OT ox) a (0 


l 


+307 “(XG “(t) + Od (x) (8) G0) g(t) 
. 1 : 1’) 


(forsee 


Now equation A.9 can be rewritten as 


i 


~ 





SG, iP (x)dm + M 9.- *(e) |] (eqn A.12) 
1 


to 


{- 
>» 


C 
=¥ a(t) a(t J 90) 9X) dm +M (6) o,(0) | 
) 


a 





: 4 
Y 4,7(t) { jy ors) dm + M 9.*(0)] 


1 


tr 


+ 


Ra Syice 4 
TG (0) 4,(0 [J O09 (x) dm +M 90 (0) | 
1"] 


-+- 


ee a 
=e LJ) Rs (x)dm + M R,-(f) | 


+ OF aK 


1 


(Ef, Rx) 04 x)dm + M R,(€) @.(£) ] 


qj 
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(for 1#}) 
From orthogonality relationship 


e 
J 01) 9x) dm + M g(x @(Q = 0 (for i*j) (eqn A.13) 


=M, (for i=j) 


where M. is the i-h mode generalized mass, then the kinetic energy of the system can 
be simplified as follows 








: 
T = “¥ q,7(t) M, + BiG, “(t) M, | (eqn A.14) 
1 = 1 





21 FR, dm + M R20) 


aD a (OFF Rl) GO) dm + MRO) OO] 


1 


From the equation 2.25, potential energy 1s 


l 
Sn 2 o- Mi: g,7(t) 
1 


(eqn Acl5) 


d. Lagrange’s equations 
Substitute equations A.I4 and A.1I5 into equations 2.19 and 
Lagangie’s equation, 
then 


2.20 to apply 





oF 


> wea ONg 20) M, + Of R200) dm + MR? (£) ] 


a uO R(x) @.(x) dm + M R,(£) 9,(0) J 
ale 9) 


Ze 


t 
= o[y q.7(t) M. + J Re) sete NLS oe | 
} 


+ Yq, (t) { FRY) af) dm + MRO 90 
1 


d BOE oe 4 > »] 
ewes flea t)M. + | R(x) VieRee (ee 
ar GIT PTE ROM # SoBe CO Amt MRO 
+287 a(t) q(t) M, 
a 
ea ee 
+ Yao (J R00 900 dm + MRO) 96) | 
l 
GU ou d aU 











oT ) 

- Sedat) vie 

Gq, 

oT a eee | 

at re ea q, (t) Mi, + O[ | R(X) p.(X) abgovens Ni RO (p( £) | 
Od, QO - 3 

d OT 


Ts eS 3a, ] = q,(t)M, + Biy R(x) (x) dm + M RYE) 9(8 ] 
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[07 


10 
Od, 





B 
= 0°, M, 9,(0) 


From all these quantities, equation A.1 becomes 


ae i ; 
O(N a(t) M, + f R-(x)dm + MR A(E) +1] (eqn A.16) 
j QO ° ZZ 


@ [a £ 
26 Yq (tv g(t) M. + ¥ a(t) | J Rd Ox) dm + M-R,(O) (6) ] = 
1 1 | 


and equation A.2 becomes 


oe 4 a 
6 | Nh R(x) (x) dm + M R,(£) @,(€) ] + g,(t) My, (eqn A.17) 
— 07 4,(t)M, + 2M, a(t) = — 250, M, 4,(v) 

fer 1,2,3,..... ,n) : 


In the equation A.16 and A.17, if we select m numbers degree of freedom of 
the system, that is, if we select » numbers of mode shape, we could have 7+ 1 numbers 


of equations. 
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DOUBLE LINK BOOM SYSTEM IN PLANAR MOTION 


a. Geometry of the system 


lo 





R(X,t) 





Figure 4.2. Parameters of the double link boom system in planar motion. 


up mass 

mass of electronics box 
length of lower boom 
length of upper boom 
angle between two links 


applied torque 


eee es 


angular displacement 
® 
Q(t): angular velocity 


R,(X) : position vector of the point on the boom in the local x-direction 
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R(Xx,{) : position vector cf the point on the boom after deformation 
W/(x,t) : deformation vector of boom 

i: unit vector of local x-direction 

unit vector of local y-direction 

unit vector of local z-direction 

unit vector of global X-direction 


unit vector of global Y-direction 


ore oa ov hou ah 


unit vector of globa! Z-direction 


b. position and velocity 
In this double link system, we have to consider the extension deformation as 
well as the bending deformation. 
The position vector of a point on the boom was expressed as a sum of R,{x) 
and W(x,t). 


Then 

Ryeet) = R(x} + W(x,t) (eqne2.9) 
where 

R,(X) mets) R(x) (eqn A.18) 


mes (Xx) 2+ RX) J 
From equation A.4+ and A.18 the position vector of the point on the boom is 


R(xt) = RADI + RX E+ VL o(x) it 20x) i] a0 (eqn A.19) 
1 


Now by differentiating the equation A.19, we obtain the velocity of the point on the 


boom. 


R(x.t) = R(x) i+ R(x) j + Y [M(x i+ OF J} Q(t) (eqn A.20) 
i 


+ Veit oF] g(r 
l 
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Substitute the tume derivative of unit vecter to the equation A.20 and simplify 


then 


R(x.t) = — OCR) + Lo 20 aly) + Foy Gv] (eqn A.21) 
1 1 
+ [R(x 0 + 8 Tos) a(t) + T oF ATi 


| i 


c. Kinetic energy and potential energy 


From the equation A.21 


R(xt) = [= OUR x) SY OMG) eo, ee q.(t) ji (eqn A.22 
i i 


+ [RY(x) 0 + OY (x) g(t) + OSC) 4,0) 13 
I l 


The dor producc of, R(x,t) 1S 
R(x.t) ¢ R(x,t) = [ — Of Rx) T S p."(X) q(t) } (eqn A.23) 
1 


we » @p.*(x) q.(t) 7 + | R(x) 8+ @ Y9;*(x) q.(t) ate Ss p."(X) q.(t) IF 


ri. ») oe 
= 9° | R,() + 2 Rx) Y oF oa(t) + Cae Ga ae 
1 1 l 


— 20f( Rw) + FY ox a0) Fox) a@] +S 1S oa ie 


1 1 l 


+ 67 [ R.7(x) + 2 R(x) V8) a(t) + OX O8(x) g,(t) a 
1 1 


20 (Ro) + Yo af} Y o(~ a(01 + (LY O20 at F 


F 
1 1 


=e 


to 


[ FRx) + RAs) } + 2E a(t) (RWW OAs) + RA) SEX) } 
1 


Se | ey d 
+ td ofp a(t) io ~ CN ox) a(t) }o] 
1 1 
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+ 26) ait) a(t) [ (R(X) + O%(x)} OP(x) — f R(x) + F(X) } OA(X) J 
l 


+1¢¥ 920) g(0}? + (LY o200 a(t) 7] 
I 1 


From the geometry of the system 
2 2 2 
R, + Ry) = Ry“(x) (eqn A.24) 


2 20) = R2 
R,7(0) + Ri-(0) = Ryn() 


Substituting equation A.24 into equation A.23 and apply the equation A.9 


then kinetic energy is 


l 


B 


T= 





ve e a) 4 + 
Q- | i Ry*(x) dm + M Ry7(€) + j,i (Y g(x) q(t)" (eqn A.25) 


1 


+ (YL oF(s) at) )*) dm + M (( L9X(0) a(t)? 
l l 


+ (Y o5(€) a(t) )*} 
1 


f | 2 
+ 2 Yat) (JC RCO O%) + Ryd F(x) ) dm 
1 


+ M (RO OE) + RE) O(8) } J 


C e ; - 
LJ ECE of) G0? + CL @700 440°} dm 
1 1 


+ 
| — 


+M{(¥ oO) g(t) > + ( ¥ P(E) g(t) Py} 
l 1 


° p . 

+ OV GOL Jt R(x) OFX) — RO) OM) j dm 
1 

+ MER (OH O70) — RYO @,(0) } 


Se c x Cie ee 
= Lal C$ (OFC) OF) — OFX) OF(0) dm 
J 
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+ M(OX(E oF(0) — 970) (0 )} | 


l 2 
+——I[. @ 
2 Ty 


Now let’s apply the orthogonality relationship 


C 
i) ot GE(X) BAX) + OFX) OF(X) } dm (eqn A.26) 


= e(H oO H+ o7(0 o*( } = (for i4j) 


=i) ( (Ona 


to the equation A.26 then 





een | if, Ry-(x) dm + M R,2(e) + ¥ g(t) M, (eqn A.27) 


a 7 
= 2 Val) (f CR Od 0°00 + RC0 9709) ) dm 
l 


~2M( RDO + RQ O2(H )} + I] 


Tez 





@ ‘s f 
~ OY GOLF (ROO) OF(X) — Ry(x) F(x) } dm 
l 
Sr Veet) eG eee Rif) (8) } 
. f 
+ Ya) FC O80 oF0x) — 9 F(x) OME) ) dm 
J 


+ M (9; (DH ae( — 97(0) @ (0) ) 3 


From the equation 2.25, potential energy is 
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U = — Yo? M, 47() (eqn A.15) 


2 1 


d. Lagrange’s equation 


Apply Lagrange’s equations A.1 and A.2 


then 
6T 
—— = 0 
0 
éT ae ‘ 
s ab eRe) +) qis(c) Vir 
35 81 | Ro-(x) dm ee 
2E ao (J CR CO O F(X) + RY(X) F(x) ) dim 
1 
+2M( ROE) + RAO OF O))+ hy 3 
e | 
ae (O19, RLCO OFEY) — RCs) O00) } dm 
a IEE SS 9 O70) — R, (£) .*(f) } 
ye ee 
72, GUUS) J 9% (xX) @- *(x) — 2 *(x) p."(8) ) dm 
j : 
+ M(9(8) eP(0) — 97(€) 9*(6) ) } | 
d OT oe ¢ 5 4 5 
—[—--] = 2(x M R,7() + ¥a.2(t) M, 
mi | aia — 
+2) 4 (of 0 OX) + Rix) O2(x) ). dm 


meee Vi ( (HOO + Ry (€) 9.” OD ae s! 


+ 20 » 4 3,(t) q(t) M. 
1 


+ 203 Gf J (Rx) Ox) + RyOx) @2(x) ) dm 
i 
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ee v , 
+ Lat) Cf) { Rd @F00 — Roo (x) } dm 
i 
+ M{R(O) OC) — RAO ONO) } 
e v : x ‘ 
+ Lait) £F, 950 0760 — 9509 90 ) dm 


+ M (OO OF — OF(0 OE} I 











éuU éU d ee 
c8 c8 dt cv 
OT wy at \ 
—— = 8 [a.(t) M, + OC) CRO) OQ) + Rx) OFC) dm 
oq, gO Wis 
~MiR(D OO + RO OH )} + 1 | 
: ZZ 
oe gv e Le eee OA, Xi) 
— q(t) lak 0, (X) D."(X) D(X) : (x) dm 
} 
~M(9,*(9 oF(0 -— o7(0 9.°(8 } J 
OT ‘ont 
—=6/(! (2 0°) + Rk Wocend 
a EF oC RG) 960 + RCO O,500 ) dm 
+ M( RO O70 - RYO) "(8 J 
3 C x y v x 
+ Veal) CF (O50 OF) — OF(0) (0 dm 
+ M (@.(0) @(0 — o2(0) @(8) } 1 + ap(t) My 
d OT oe C ; | a 
ene re =| i CRO) OO) + RCS) 708) ) dim 


~ M( RO) (0) - RO O, (0) } 
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£ ie : 
+ . q,(t) { Je P(X) O(X) — OF(X) O, (0) dm 


+ M (OD OF(D — O2(9 O(9 } 1 + aC My 


OY 4, (CJ, ‘ @.x) Oy (X) — 9(X) O,%(x) } dm 


l 


+ MLO (OOO — 927( 0,0] + a,(0 My 





3 
= O°), My (0 





Plug these quantities in the equation A.1 
then 


a ; 
0 | J, Ry(x) dm + M R,*(€) + 2 g.7(t) M. (eqn A.28) 
1 


wl | 
+ 2V a(t) CJ CRO) OF) + RCo) @F(x) ) dm 
i 


+ 2M( R (C) .*(C) ES RO p.>(£) at an! 


+ 26 eq 3.(t) q.(t) M, 


1 


at J (R(X) O A(X) + R(x) O2(X) ) dm 


+ 
tI 
se) 
bss ie 


*2M( RAD OH + RO OO )} + Tp 0°} 


+ 


_ we ae | 
LGD £6 OFC) FO) — OFC) GAH } dm 
J 


M(OA(9 oD — 92(8 90 | = 


+ 
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and the equation A.2 becomes 


eet | 
6 [ i CRiolroney: R(x) @,*(x) ) dm (eqn A.29) 
+ M( RO 9,7(8 - RAO 9, (©) } 


TOY SOG; (8) 0 OO a aes: 


+ 
ee id 


~ M £@*(0) (0) — @5(0 9:(0) } 1 + (0 My 


: : : | 
Donley Mee na ROD OO) + RCO O,"(X9) ) dm 


+ MER (OOO + RO O5(0 D3 | 


ti 


Dias hae e a8 ) (x) > “(x) ( Miy\ ig 
2 = HOU (X} P(X) — @F(X) Oy *(x) } dm 


Mi e(0) 9, (€) = 0p," (€) oO (€) Piles oO”, mn q(t) See . Ww Mi qn) 


e. Results : 

Table 5 and 6 show the eigenvalues of Aluminum Alloy and Composite 
material boom from. NASTRAN simulation results and Figure A.3 shows first two 
mode shape of the double link boom in planar motion. As shown in Figure A.4, 
applied torque to maintain 15 r.p.m of rotating speed increased to one hundred times 
of the torque that applied to the same boom in 3-dimensional motion because its mass 
center 1s far awav from rotating center. Angular displacement change with time was 
almost same as 3-dimensional case. Displacement in x and v-direction was much bigger 
than that of 3-dimensional motion because tip position is far awav from the rotating 
center and concentrated mass attached to the tip. 

While rotating speed increases. boom deflects negative v-direction and positive 
x-direction until reaching its equilibrium position. After cutting off torque. boom 


remains equilibrium position and oscillates with some amplitude. Magnitude of 
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deflection at equilibrium position was about 6.2 inchs and slope at the tip position was 
2.3 °. These values are much bigger than that of 3-dimensional mction but it is quite 
natural because of its big radius of rotation. 

Figures A.5 - A.12 show angular displacement, angular velocity, generalized 
displacement, deflection in each direction, magnitude of deflection and slope at the tip 


position respectively. 


ASGIEE 5 
Ineo EIGEN VaLUES OF ALUMINUM ALLOY (2D) 


= PO OSONGl > 00, 5.67997 =O] 1 .QOOQO0E + 00 P27 205s 0 
Bao 55 + Ol Ts oI DOI UU 1 .QOOQOOE + 00 D.22ehok 702 


TABLE 6 
Re etal VALE ES OF COMPOSITE MATERIAL ( 2D } 


mode hors evcies eae generalized 
no. W); mass( sulfness( Kk.) 


4. | 4.461328E+00 | 7.100423E — 01 1.QOQOO0E + 00 1.990344E + 01 


2.246115E+0901 3.574804E + 00 1 .QOOOOVE + 00 5 US5055 5 7, U2 









generalized 
suffness(k.) 


generalized 
mass( M-) 
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Figure A.3 First and second mode shape. 
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Figure A.J. Applied torque vs. ume. 
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Figure A.) Angular displacement vs. time. 
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Figure A.6 Angular velocity vs. time. 
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First mode generalized displacement vs. ume. 
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Figure A.8 Second mode generalized displacement vs. time. 
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Figure A.9 Displacement in x-direction vs. ume. 
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Figure A.10 Displacement in v-direction vs. time. 
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Figure A.11 
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Magnitude of deflection at tip position vs. tune. 
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APPENDIX B 


DETAILED DERIVATION OF LAGRANGE’ E UNOS FOR THREE 
DIMENSIONAL MOTIO 


Applv Lagrange’s equations 2.7 and 2.8 into equations 2.28 and 2.25 





then 
La 
G8 
6T . ec 
ae = PLL a M, + J R(x) dm + M R,(E) 
ae 
+ 2h at) (J, RO (x)dm + M R,(£) p.*(€) } 
i | 
ast 
-2tf (Yo, det yo? + ky Mi 
l l 
ma () a(t) ( Fax) C8) dm + M @;*(€) 97(E) } 
ry 
Sea (oa 2 RA (x) O- h(x eins M R,(€) 9 (£) } 
d GT “a a eas m, 
—— {—y—] = 0[3S a.~(t) M. + {| Ro“(x) dm + MRE 
pa era [ S a,7(t) M; iF (x) dm he 


l 


2 ad (0 (J RAW) 900 dm + MRO (6 } 


Je 
= 2¢) (> OAx) a (0)? + POD Oey eee 
C1 ; (ld 


ra 2054 (t) q(t) M, + j RCo .*(x) dm + M R,(0) (8) } 
1 


Le 
Y 4) CF 091%) Ox) dm - M 9%) O18 } | 


e@ e 
DY at (L a,(t) { \, @(x) OF(x) dm + M (0) @ "(8 
l 
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Sar | 
— J0F 0) 9/809 dm — M O70 ONO 5 


C _ ; 
- FR, 9) dm = MRO 950] 


OT : 
= @[ qit) M, + f R(x) o(x) dm + MR,(O) 9°(0) 
04, { ] 0 : 1 ax 1 
c Vs 
— 26 a(t) £f 9700 O° dm + M 9,70 9,7(8) 3] 
~ +} gi! { } 
1 
+9vc et ‘(x)dm + M@,*(0) 5(€) 
— q(t) | gen (X) Q(X) OO Nee he oF (C ta | 
l 
C Vv 5 a 
- J}: i) @,*(X) dm - M 9.*(0) o(0 
4 (t) M 
— rns 
6G, t h 
0S iq Te, dm + M o.%(0) o (0 
a q(t) Ty (X) O° (X) mM : ?. ( ) Oi" ) 
l 
e Vv x ay x 
= RO (X) P, (<jdm- i p.(C) He (€) 
> a 
+ O[ \, R(x) @,"(x) dm + M R(€) p90) | 
d | OT , ie 
emis. . | - Gt). 
dt 6q, Wnt) Ma 


ee t 
+ OLX a(t) (J OF) O80) dm + MONO 9,0 


1 
y x iY x 
= } 9 @, (x) dm - M @7(0) 9, "(0 } 


iP 


C | , 

a pRalX) On'(X) dm + MRO O,°(8) | 

Oy Wa) M @X(£) o:5(E 

+ OD ad LJ OR) ys) dm + MONO OA 
] 


or | ; | 
= jr) o,*(x)dm- M 92(0 9, %(0 | 








a 
c8 66 dt G6 
ca Bie allo 
= (05 t 
6a, h h Uys 


-_ 


Now plug all these quantities in equations 2.7 and 2.8 
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APPENDIX C 
NASTRAN PROGRAM FOR DYNAMIC ( MODAL ) ANALYSIS 


1 DYNAMIC ANALYSIS IN PLANAR MOTION 


S 

§$ Double link flexible boom for dvaamic analvsi tia 14 grid points ) 
S in 2 dimension space with tip mass (37.5 16. ) linkI714 ft,link2: hy ft 
S  Material:aluminum alloy 

ar lctinses kaeek dun sansncaaase 

S 

meee ECU LIVE CONTROL DECK 

errr errr rere ree 

S 

id kang,dynamics SN-ROSS 

sol 3 

time 10 

diag 8,13 

ne 


tag se ESS HERES Ee NCEE 6 2s ENE 2 OE AE A Oe Nes 


S 
: ese CONTROL DECK 


Ste 20 20s 8s 25 SC 6 25 IC He OR Des ae 2 Ae ana ae OO oe 


title= Modal analysis of double link flexible 
title= boom in planar motion( Aluminum alloy) 
echo = both 
os 120 

c= 101 


ae all 

put(plot) 

: Biter Hestran 

axes Z.X.V 

cscale = {1.8 

view 0..0.,0. 

as ee 14.0 by 10.0 


find scale, origin l.set | 
plot modal deformation _0,set L,origin l,shape 
maximum deformation 5 


ean bulk 

Ipeiaseai cts cee 020 oe nee che 0 EC OR 2E He 22 fe 2G he 2 abe ae She ae aft he 21 ale 2 she a oie 
S 

| BULK DATA DECK 

S** NG mesiee ah Ope arm O40 2 ME She eee che ee ogc Sie ook cae wie SHE ace See Saltese cee ote He ote SK OK 
S 

: Define new coordinate system for convinience 
peord2c,!,,0.,0.,0.,0.,0.,1., + 23 

emia 5, 1.,0., 0 

S grid(node) data 

grid,1,,0.,0.,0 

= *(1),= ,°(24.),= = 
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Hu ti 
- 
Ch 
iC 
— 
li co 
mts 
& 
CN 
ore) 


40° 5. 


lement data 


bar, ar,21.102,1,2,0..0.,1. 
1),=,"(1),"(1),= 


rt 
Change these boom os when material 
changes AL. to COMP 


S --- DATA FOR ALUMINUM ALLOY --- 


QO UFUNA 
— cb 


spcadd. 
eigr, 1 20,mgiv,, 
Sparam, aUutospc,} yes 


S --- DATA FOR COMPOSITE --- 


S 
Pear iO2 105-2. OlGse ae lls oles 3.8902 
mat!.103, 05) 5562: 
.103.14..9 176-2 
6 


Cok a 
SOCinu 2 2. thru.14 
spcadd. | ee 
e1OT: | 20 M1e1V 
Sparam,autospc, ves 
enddata 


Art 


I 
} 
3 


reer 
cy 


ie 
S-4 
Ol, 


aay OT Oey en 
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2, DYNAMIC ANALYSIS IN 3 DIMENSION SPACE 


Sy 


mpd link flexible boom for dynamic ane alvsis( 14 grid points ) 
3 dimension space with tip mass (37.5 Ib: ) 
HAR I: eae ft tink 2 
angle between LI and x-axes is 70,L1 and L2 is 126 degree 
Material: Aluminum alloy 


te Me os fe ste ote ole ne ois ate fe ate aie ais ote fe oft nis fe 20 ae ale ole ae ofs st ote cis ois als ole 2s ote fe oie 


Pec @leiIvE CONTROL DECK 


ste afc ofe ate ale fe ate ate aie ae ae a ale abe ole ae fe ate fs a a aft ale Xe ale aoe fe ale ois as 


L, 
sf 
% 
x 

vi 


Se el git om 


id kang,dynamics Snross 
sol 3 

time 10 

diag 8.13 


cend | 
Se ee Me Ba ee aGae a ea pe ole HE ole ote ato he ots ote ste ste cle is Xe ahs fe ote ote ofe safe of 


@yse CONTROL DECK 


’, 2. 4, ?, o U f, 
ie of os 3s Me ote sis oie ots als nie ois ois ale ofc oie ote ole oie ae ate of oie ols rie ote ots oft ois ole He ok 2k ok oie 


Ae. 


title= 2 link flexible boom in 3 dimension(Al. COM) 
echo = both 


isp = all 
output( plot) 
plotter nastran 
cscale = |. 
view ()..0..0. 
paper size 14.0 by 10.0 
set l=all 
find scale. origin l.set | 
maximum deformation 10 
aes my .X,Z 
plot modal deformation 0,set l,origin l,shape 
axes X.Y.Z 
lot modal deformation 0.set l.origin l.shape 
egin bulk | 


se Xe ole fo ate ois ale ole ole ote Me Ne ste ofs ate ole ale slo ste fe ote aie eo als 


BOLK DATA DECK 
308 HS HE SS NS HS os Ss He aS os he os oN ae oie 2s I Ne ats ole aft ate cle Nt ats oie ole fs ale cic ie 
Define new coordinate system for convinience 


Bene, 1.,0.,0.,0..0..0.,1., + 25 
peer ),U.,.906/ 


la a 


5 
grid 0.0.9 7 
1). = ,*(8.2085), = ,*22,5526 


<a 


S 
S element data 


[25 


QO UVMA 
o G 
2) 


ir 


a 


il 
Change followin aes when the material 
changes from A 


S --- DATA FOR ALUMINUM ALLOY --- 


pbar,102,103.1.2101.1.2447,1.2447,2.4894 
matl.103.1.01 > pe ae 2 .3362-4 
conm2. 103. |- 2 
spcl, ese 
ies: 
101, 
elete) 20: mgiv,, oe 
Sparam, autospe. yes 


: --- DATA FOR COMPOSITE --- 
aban Whe aCe ne af wt l. 2447. 1.2447 ,2.4894 
2.5302-4 


(SAV IAY tl tl 


ee 


Pants 


eee ‘autos aaa yes 
enddata 
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APPENDIX D 
DSL PROGRAM SOLVING THE DYNAMIC EQUATIONS OF MOTION 


1. PROGRAM OF DOUBLE LINK FLEXIBLE BOOM IN PLANAR MOTION 


i This program solves the dynamics of a dobule link flexible boom system 
> in planar motion. The applied torque ( t ) was changed to see the effects 
= of the torque and also damping coefficient ( ¢ ) was varied 0.0 to 0.5 %. 
. This program automatically calculates the slope and deflection of the tip 


t position in each direction ( local x, y, z ) from simulation results. 


; N 
FL1; length of lower boom 


- number of discritized boom element 


* —_ FL2; length of uppper boom 

- vi: tip mass 

* RHO ; mass per unit length 

* DX; differential length between each grid point 
RX ; local x distance from origin 

* TAO; applied torque 

= OMGI ; Ist mode natural frequency 

* OMGz2 ; 2nd mode natural frequency 

* PX1 


PY1; Ist mode shape-displacement in local v direction 


Ist mode shape-displacement in local x direction 


wee 


- PZ1.: Ist mode shape-displacement in local z direction 
7 PX2; 2nd mode shape-displacement in local x direction 
, PY2; 2nd mode shape-displacement in local v direction 


“ PZ2; 2nd mode shape-displacement in local z direction 
ALP; angle between upper boom and local x-axes 


RXIL ; Ist mode shape-rotation at the tip position W.r.t x-axes 


ey 


oe RZIL ; Ist mode shape-rotaticn at the tip position w.r.t Z-axes 
RY2L ; 2nd mode shape-rotation at the tip position w.r.t y-axes 
‘9 RZ2L ; 2nd mode shape-rotation at the tip position 

: TH; angular displacement 

¥ THD ; angular velocity 


: Ql; Ist mode generalized coordinatest 


is QID ; time derivative of lst mode generalized coordinates 
5 2; 2nd mode generalized coordinates 
a Q2D ; time derivative of 2nd mode generalized coordinates 


* ZETA ; damping coefficient 

= RIRZZ ;mass moment of inertia for R.F electronic box 
RIBZZ ;mass moment of inertia for boom and tip n mass 
* Cl; starting up time coefficient 

- C2; magnitude of applied torque 


i C3; ‘required time to reach requinedakel ov 


i 

Me ae she HHS she ate ote Ns ake ots fe ahs ote ahs ots fs ote ote hs He ots afc oe oe ofc fs nis ake ote ots ts ote ots of ote ofe ots oie ote ofe sts ake ofe ote fs ts ate fe fz ote oft ote ots oe mie ofs ote ote oe ote ote ote ofs ote ste ote ate ats ots ahs ake aks ofs ahs ote ote ots ahs oe ots 
we 

* SIMULATIOMS OF DOUBLE LINK FLEXIBLE BOOM IN PLANAR MO Gitges 


Be He oe Me ee oe Ms ote oie ae aks of ate He ois ote aie ote ois Me aie oe ole ois ate ote ote ofa aks cis ole ote ote ake sie os ake Me he ote fe fe fe ale ols ots ole ole oe ahs oe ate ols ake ote oe Ne ris sfe oie ote ake afe oft oe ake oft Me cis ot oie ie Ne Ee oe 2s oe oe 
te 


TITCTE SOLUTION OF SIMULTANEOUS DIFFERENTIAL SE OG aes 
TITLE FOR DOUBLE LINK FLEXIBLE BOOM IN PLANAK iO aie 


FIXED IER. IPVT.N. 1 7 
CONST FLI= 168..FL2= 144..RHO = 3.0690E-04.N = 13.DELT=.04 
PARAM Cl1=40. (C2= 1000. C3= 19.65 
PSA 
TM =9.7176E-02 
ALP =34.  _. 
RIRZZ = 31.0559 
RAG = bop rise on 
; RZ?IL =-3.149482E-02 
* BY CHANGING THE ZETA VALUE THE DAMPING FORCE 
* WILE BE CHANGED 
© ZETA=0. 
ZETA=0.002 __ 
. ZETA =0.005 


Ie 
ne 
on 
y 
oY 

O 
2 

& 


ses @@E ee ee @ PF 2 SF SSF FS SFT FFF FFFFW Bt SF STF FT FIsF SFT FSF SF FFs Fs FH FTFFFT FI BGFFZsF BFF BFF SFT FT FFF FF FSO et Fee eT OeneZ 2 e222 ee e22 @ee @ 
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